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ABSTRACT 



A popular model for the circumstellar disks of Be stars is that of a geomet- 
rically thin disk with a density in the equatorial plane that drops as a power 
law of distance from the star. It is usually assumed that the vertical structure 
of such a disk (in the direction parallel to the stellar rotation axis) is governed 
by the hydrostatic equilibrium set by the vertical component of the star's gravi- 
tational acceleration. Previous radiative equilibrium models for such disks have 
usually been computed assuming a fixed density structure. This introduces an 
inconsistency as the gas density is not allowed to respond to temperature changes 
and the resultant disk model is not i n vertical, hydrostati c equilibrium. In this 
work, we modify the bedisk code of ISigut fc Jonesl (120071 ) so that it enforces a 
hydrostatic equilibrium consistent with the temperature solution. We compare 
the disk densities, temperatures, Ha line profiles, and near-IR excesses predicted 
by such models with those computed from models with a fixed density structure. 
We find that the fixed models can differ substantially from the consistent hydro- 
static models when the disk density is high enough that the circumstellar disk 
develops a cool (T < 10, 000 K) equatorial region close to the parent star. Based 
on these new hydrostatic disks, we also predict an approximate relation between 
the (global) density-averaged disk temperature and the T ef j of the central star, 
covering the full range of central Be star spectral types. 



Subject headings: stars: circumstellar matter - stars: emission line, Be 



-2- 



Introduction 



Be stars are non-supergiant B stars that currently show, or have shown in the past , 
emission in one or more of the hydrogen Balmer lines, typically Ha (IPorter fc Riyinusl 120031) . 



Other common characteristics of Be stars include an infrared excess (Waters 



1986) and linear 



conti nuum polarization at the level of approximately one percent (IWaters fc Marlborough 
19921 ). These observations can be explained by the pres ence of a non-sphe rical distribution of 



circumstellar material surrounding the central B star (IWood et al.lll997). Be stars have re 



cently been resolved interferometrically in the optical (jQuirrenbach et al. 



19931 ; iTycner et al. 



20051 . 120061 ) to conclusively show that this cir cumstellar material is in the form of a thin equa- 
torial disk, a form originally championed by lPoeckert fc Marlborough! (119781 ). Typically the 
spectra of Be stars are consistent with this thin disk in near Keplerian rotation in which the 
equatorial density drops as a power-law in radius. In such models, the vertical structure of 
the gas (perpendicular to the disk) is often assumed to be in isothermal hydrostatic equi- 
librium. In this case, the overall disk density, expressed in cylindrical co-ordinates (R, Z) 
where the Z direction is parallel to the rotation axis of the star, is 



Pi /?• Z ) />„ ( ^ 



X 



Here R* is the stellar radius, p Q is the density of the inner edge of the disk in the equatorial 
plane (often termed the "base density" of the disk), and if is a vertical scale- height that 
depends upon the distance, R, the stellar mass, M*, and a disk temperature, T Q , assumed 
valid for all Z at that distance. Such a density model balances the vertical component 
of the star's gravitational acceleration (parallel to the rotation axis) with the pressure set 
by the temperature T . If T Q is a constant for all radial distances, the weakening of the 
vertical gravitational acceleration with R causes H cx i? 3//2 and the disk flares with increasing 
distance. Here p a and the power-law index, n, are the parameters that determine the overall 
density of the disk (together with the assumed value for T Q which is often simply set to a 
constant fraction of the T e g of the central star). Typically one fixes these parameters by 
comparing to the observation al constraints cited abov e. By matching t he in frared excess 



of a wide sample of Be stars, IWaters et al 



9871 ) and iDoughertv et a 



19941 ) find power- 



law indexes in the range of n = 2.0 to 5.Qj. Theoretic ally. iPorterl (119991) no t es tha t an 
isothermal, viscous disk requires n > 3.5 for outflow, and lJones. Sigut fc Porter! (120071 ) find 
that a hydrodynamical simulation of outflowing viscous disks in which the thermal structure 



1 Both of these analysis are based on the disk model of Waters! (1986) which assumes a density distribution 
of the form p(r) = p (R^/r) n in a disk of finite opening angle. Here r is the distance from the centre of the 
star. This density distribution thus differs from Eq. [TJ 
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of the disk is taken into account (following I Jones. Sigut &: Marlborough! 120041 ) predicts n in 
the range of 3.0 to 3.5 in the inner portion of the disk, R < 10 i?*, where hydrogen spectrum 
likely forms. 

While empirically probing the temperature structure within Be star disks is a difficult 
observatio nal challenge, several theoretical models of the te mperature structure of such disks 
now exist (jCarciofi fc Bjorkmanll2006l ; ISigut fc Jonesll2007l ). Typically, for low densities, the 
disks are nearly isothermal. However, as the density rises (i.e. as p Q is increased), a cool 
equatorial region develops close to the star and the disk gas becomes far from isothermal 
in the vertical direction, even at a single R. Such models are fundamentally inconsistent 
because the derived temperature structure can no longer be consistent with the assumption 
of a vertically isothermal gas as required by Eq. [TJ As observational diagnostics can be 
very sensitive to the gas density in the disk, such an inconsistency may have observational 
consequences in terms of the predicted emission line profiles and stren gths, the predicted IR 
excess , and the predicted linear polarization sign ature. F o r exa mple, ICarciofi fc Bjorkman 
(120061 ) suggest that the fundamental limitation to IWaterd (119861 ) technique for determining 
disk density from the slope of the IR continuum is the a priori assumption of disk geometry 
via the adoption of a constant disk opening angle. 

It is the purpose of the present work to eliminate this inconsistency between the tem- 
perature and density structure of the disk. Radiative equilibrium models are constructed in 
which the vertical disk density structure is in a hydrostatic equilibrium consistent with the 
computed temperature distribution. 



2. Theory 

We assume an axisymmetric disk described by the cylindrical co-ordinates (R,Z). We 
shall also assume that the disk density in the equatorial plane (Z = 0) is a known function 
of R of the form 

p{R, 0) = p, 

where R* is the stellar radius; p Q and the power-law index n are free parameters that fix the 
density structure of the disk. Thus at any location R, the density p(R, 0) is known, and we 
shall denote this density simply as p(0). If the gas at this location is in vertical hydrostatic 
equilibrium, then the vertical pressure gradient must satisfy the equation of hydrostatic 
equilibrium, 

dP 
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Here g z is the vertical (or z- component) of the star's gravitational acceleration at location 
R, namely 

9z = ° M * (R 2 + Z 2 ) 3 / 2 ' (4) 

where M* is the mass of the central star. In all cases, p Q in Eq. [2] is so small that the 
mass of the disk is completely negligible compared to that of the star. The pressure can be 
eliminated from the hydrostatic equation as the gas is assumed to obey the perfect gas law, 

P kT . (5) 



pm H 



This introduces the vertical temperature distribution into the problem. In this equation, P, 
p, T, and p are functions of Z (the mean-molecular weight, p, depends on Z because it is 
determined by the ionization state of the gas). Using the perfect gas law to eliminate the 
pressure from Eq. [3J and using Eq. HI we find that 



where ot(Z) is the function 



- P tz = ~ a(z) (B TzW ' (6) 

Using the boundary condition that the density at Z = is equal to p(0), Eq. [5] can be 
numerically integrated assuming that the functions T(Z) and p{Z) are known. Typically a 
numerical solution is required as the temperature T(Z) is available only at a fixed number 
of vertical grid points v ia a radiative equilib rium solution. In the sections to follow, we shall 



use the bedisk code of ISigut fc Jonesl (120071 ) to compute the thermal structure of several Be 
disks so that the vertical hydrostatic equation can be numerically integrated at each distance 
from the star. 

The analytic solution of Eq. [61 seen in Eq. [H is obtained when the vertical temperature 
variation can be replaced by a constant temperature T a (with the mean-molecular weight 
also assumed to be a constant, p Q ). In this case, Eq. [7] is just the parameter a , and Eq. [H] 
is readily integrated to give 



P (z) = p(o)e y R 

In the case of a geometrically thin disk in which Z/R <C 1, applicable to the Be stars, this 
result simplifies to 

p(Z) = p(0)e"^(^) 2 . (9) 
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Comparing to Eq. [TJ the scale-height H is given by 



H=J^, (10) 

which gives the cited H oc R 3 ^ 2 disk flaring with distance. Thus the assumption of vertically 
isothermal gas in the Be disk is required to reproduce the density structure of Eq. [TJ 



3. Computational Procedure 



To obtain the tem p eratu re structure of various Be star disk models, we use the bedisk 
code of ISigut Jonesl (120071 ). This code has been succes sfully used to interpret a wide 



range of B e star observational signatures (for example, see iTycner et al.l 12008 



Jones et al. 



20081 . l2009l l The bedisk code solves the radiative equilibrium problem for the disk gas by 
balancing the heating and cooling rates computed for a user-specified set of atomic models. 
Energy input to the disk is assumed to come from the radiation of the central star. This 
photoionizing radiation field is represented as a sum of the direct component from the central 
star and the diff use component fr o m the disk itself , which is treated by the on-the-spot (OTS) 
approximation fjOster brockill989l ) . ISigut fc Jonesl (120071 ) present results for a direct (but still 
approximate) treatment of the diffuse component to demonstrate its affect on the thermal 
structure of the disk. They find that even for a dense disk with p Q = 5 • 10~ n gem" 3 and 
an R~ 2 5 radial drop-off, the additional heating provided by the diffuse component increases 
the temperature in the inner equatorial plane by only 10-20% and conclude that the OTS 
approximation gives good results. 

Limb darkening of the central star is accounted for in the direct photoionizing radiation 
field, but the gravitational darkening and geometric distortion implied by rapid rotation is 
yet to be implemented. Heating of the gas is via photoionization and collisional excitation. 
Cooling of the gas is via recombination and collisional de-excitation. The atomic level 
populations required to compute the heating and cooling rates were found by solving the 
requisite statistical equilibrium equations. The transfer o f line r adiation is handled through 
the escape probability approximation. See ISigut fc Jonesl (120071 ) for further details. 



To achieve a solution that is both in radiative equilibrium and in a vertical hydrostatic 
equilibrium consistent with the temperature structure, some modification to the BEDISK code 
is required. The initial density distribution is simply that of Eq. [TJ bedisk then obtains 
the disk temperatures on a fixed grid specified by the cylindrical coordinates (R4, Zij), with 



1,2, 



n, 



and j = 1, 2, 



n. 



where R4+1 > Ri and > Z^j. The solution starts 



closest to the star (at i — 1) and statistical and radiative equilibrium is solved for each Z\ 
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for j — n z , n z — 1, . . . , 1, i.e. by proceeding down towards the equatorial plane. Then Eq. [6] is 
numerically integrated using the density in the equatorial plane as the boundary condition. 
As a result of this integration, the total gas density for i — 1 and all j = 1,2, . ..,n z is 
updated. As the gas density at each Zij has been changed, radiative equilibrium must 
be resolved to obtain a new temperature distribution. This process is iterated until the 
maximum fractional change in the gas density drops below 1%. Once this has happened, the 
calculation proceeds to the next radial distance in the disk, i = 2, and so on. 

If, on average, N p density iterations are required to meet the convergence tolerance, 
then the entire calculation for consistent vertical hydrostatic equilibrium takes N p times 
that required for the radiative equilibrium solution with the fixed density structure of Eq. [TJ 
As N p ~ 5 — 10, a significant lengthening of the execution time of bedisk occurs. Finally, we 
have found that the various iterative estimates for the density often tend to oscillate around 
the true value, and it is generally advantageous to update the grid densities with a relaxed 
estimate that is the average of the old value and the new predicted estimate. 



4. Results 

To demonstrate typical differences between disks computed with a fixed, isothermal 
density structure (referred to as "fixed" models) and those computed with consistent verti- 
cal hydrostatic equilibrium (referred to as "hydrostatic" models), we consider pure hydrogen 
and helium envelopes surrounding the three central stars with parameters found in Table [TJ 
These parameters represent central stars with spectral types of approximately BO, B2, and 
B5. While the bedisk code can handle multiple atomic models, and thus can compute 
the radiative equilibrium solution for a gas with a solar composition, such calculations are 
much more computationally intensive than the hydrogen/helium models considered here. 
Nevertheless, the hydrogen/helium models are completely adequate to illustrate the main 
differences between the fixed and hydrostatic models over a wide range of disk parame- 
ters. For the atomic models, we have adopted the 15 level Hi and 13 level He I models of 



Sigut fc Joned (120071 ); Hll and Hell/ill were represented as single levels. While the inclusion 
of helium has only a small affect on the thermal solution, it allows for a more realistic mean 
molecular weight for the gas. The total helium number density was assumed to be 0.1 of the 
total hydrogen number density. 

For all calculations, the disk was represented by N r = 84 and N z = 50 grid points with 
Rmax/R* = 50. The power law density index for the equatorial density (Eq. |2J) was taken to 
be n = 3.5 in all cases. 
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In the discussion that follows, we first outline the detailed differences between the pre- 
dicted temperature and density distributions in the disk models. Next we consider how well 
the fixed models can predict the density-weighted disk temperature given an appropriate 
choice for the T parameter. Finally, we examine the differences in selected observational 
predictions of the models, namely in the Ha line profile and in the near-infrared excess, both 
of which have been used by previous investigators to determine the density structure of Be 
star circumstellar disks. 



4.1. Disk Temperatures and Densities 

Figure [1] compares the predicted density structure of fixed and hydrostatic disks com- 
puted for the BO model with p Q = 5 • 10~ n gem -3 . In the case of the fixed model, an 
isothermal temperature of T a = 13, 500 K was chosen for Eq. [1] (the reason for this choice 
will be discussed below). The lower-left panel of this figure shows the difference in the log- 
arithm of the predicted density, and the lower-right panel, these differences as a histogram 
that includes all of the grid points. The hydrostatic model is generally more concentrated 
towards the equatorial plane than the fixed model, and this is particularly clear from the 
histogram of differences. This central concentration is a result of the cool, equatorial region 
that develops in a such a high density disk (see next paragraph). 

Figure [2] compares the predicted vertical disk scale heights between the fixed and hy- 
drostatic disks. Scale heights for fixed disks follow from Eq. [10] and are shown for isothermal 
temperatures ranging from 8500 to 17500 K. Scale heights for the hydrostatic disk are found 
numerically by locating the point above the disk at which the density falls to 1/e of its 
value in the equatorial plane. The cool equatorial region that develops close to the star 

Table 1. Adopted parameters for the central B stars. 



Parameter 


B0 


B2 


B5 


Radius (R©) 


10.0 


7.0 


5.0 


Mass (M ) 


17.0 


12.0 


9.0 


log 10 (L/L e ) 


4.5 


3.8 


3.0 


T eff (K) 


25000 


20000 


15000 


log(#) 


4.0 


4.0 


4.0 
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(see additional discussion below) is clearly reflected in the hydrostatic disk scale height for 
R/R* < 5. Beyond this, the scale height rapidly rises, and by R/R* ~ 10, it closely matches 
the fixed, isothermal model corresponding to T Q = 13, 500 K. Thus in the region close to the 
star, R/R* < 10, the scale heights and their variation with R are not well represented by 
any of the isothermal models. 

One additional subtle point is that the hydrostatic disk is somewhat less massive than 
the fixed disk. For the fixed model, the mass follows directly from Eq. [T] and the adopted 
parameters. However, for the hydrostatic disk, the mass follows the density adopted in the 
equatorial plane (Eq. [2]) and the hydrostatic equilibrium solution. In this case, the mass of 
the fixed disk is 2.78 ■ 10 -8 solar masses whereas the mass of the hydrostatic disk is about 
20% less at 2.16 ■ 10~ 8 solar masses. 

Figure [3] compares the predicted temperature structure of both disks. As shown in the 
lower panels, there are significant differences in temperature. Due to the large initial density 
(p = 5 ■ 10 -11 gem" 3 ), both the fixed and hydrostatic models have a cool, equatorial zone 
close to the star. Above and below this cool zone are hotter sheaths which can still be directly 
illuminated by at least part of the central star. Large temperature differences between the 
fixed and hydrostatic models can result from the differing locations of these hot sheathes as 
the hydrostatic model is more concentrated in density towards the equatorial plane. There 
is also a significant temperature difference between the two models at the location of the 
optically thin gas far above the equatorial plane. This is a result of the very large density 
difference between the fixed and hydrostatic models in this region as illustrated in Figure [TJ 

In this comparison, one might question the choice of T Q = 13, 500 K adopted for the 
fixed model. To address this point, an additional set of fixed calculations was performed 
by varying T Q from 9,500 to 17,500 K with all of the other model parameters held fixed. 
Figure H] summarizes these results by giving histograms of the differences in temperature 
and in the logarithm of the density over the entire grid. As can be seen from Figure HI the 
density distribution is best represented by the coolest model, T Q = 9, 500 K. This is not 
surprising as p Q is large enough that a cool equatorial region develops with temperatures as 
low as 8, 000 K in some regions. Conversely, Figure H] shows that the temperature is best 
represented by the hottest model, T Q = 17, 500 K. This result, however, is a bit deceptive. 
As will be shown in the next section, it is the coolest model which does the best job in 
reproducing the observables. The large temperature differences in the T Q = 9, 500 K model 
often result from the misalignment in the height of the cool gas and hot sheaths. 

The previous large differences between the fixed and hydrostatic models reflect the poor 
job the assumption of an isothermal T a does in representing the density structure of a disk 
that develops a cool equatorial region. Figure [5] illustrates the development of this cool 
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Fixed Hydrostatic 




2 4 G 6 10 -3-2 -1 1 

R / R„ A. log a (hydro-fixed) 



Fig. 1. — The density structure of a fixed (upper left) and hydrostatic (upper right) disk 
with p = 5 • 10~ n gcm~ 3 and n = 3.5 for the BO model of Tabled! The lower-left panel 
shows the differences in log(p) directly, while the lower-right panel gives a histogram of these 
differences over all 4200 grid points. 
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0.2 




0.02 - 



2 4 6 8 10 12 14 16 18 20 

R/R* 

Fig. 2. — Vertical disk scale height as a function of radial distance for a BO model with 
p = 5 i0~ n gcm~ 3 and n = 3.5. The solid lines are scale heights for fixed models (EqnsQ] 
and [TUl) with isothermal temperatures ranging from 8500 to 17500 K in steps of 1000 K 
(moving upward in the figure). The filled circles are the scale heights of a hydrostatic model 
of the same p Q and n found by locating the vertical height at which the density dropped to 
1/e of its value in the equatorial plane. 
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Fig. 3. — The temperature structure of a fixed (upper left) and hydrostatic (upper right) 
disk with p a = 5 • 10 gem and n = 3.5 for the BO model. The lower-left panel shows 
the temperature differences directly, while the lower-right panel gives a histogram of these 
differences over all 4200 grid points. 
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10000 



10000 



10000 



A log p (hydro-fixed) A T(K) (hydro-fixed) 

Fig. 4. — Histograms of the density (left panels) and temperature (right panels) differences 
between fixed and hydrostatic models for p Q = 5 • 1CT 11 g cm" 3 and n = 3.5 for the BO model. 
The fixed models assume T Q = 9,500 K (top row), 13,500 K (middle row), and 17,500 K 
(bottom row). 



region, for the BO spectral type, as the overall disk density is increased. Plotted is the 
density-weighted temperature, 



< T Disk (R) > 



J Zmax p(R, Z) T(R, Z) dZ 
Jo m ~p(R,Z)dZ 



(11) 



as a function of radial distance for six hydrostatic models with initial densities, p a , ranging 
from 5.0 ■ 10 -13 to 1.0 • 10 -10 g cm -3 . All models assumed n = 3.5. This figure illustrates that 
models with p a < 5.0 ■ 10 -12 gem -3 , do not have an extensive cool, equatorial zone near the 
star where T(R, Z) < 10, 000 K whereas denser models rapidly develop such a region. This 
lack of a strong temperature gradient for the less dense models suggests that the isothermal 
approximation of Eq. [TJ is adequate to represent their density structure. This is borne out 
by Figure [6] which compares the disk temperatures of a low-density (p Q = 1.0 ■ 10 -12 gem -3 ) 
hydrostatic model with a fixed model computed with the same p Q and T a = 13, 500 K. There 
is little temperature difference between the two models. Figure [7] compares (as histograms) 
the predicted disk temperatures and densities for three choices of the T Q parameter, 9, 500, 
13, 500 and 17, 500 K, with the hydrostatic model. As is clear from this figure, T Q = 13, 500 K 
does an adequate job of reproducing both the densities and temperatures in the disk. 

The general conclusion is that the isothermal, fixed T Q models are appropriate for low- 
density disks provided a reasonable choice for T a is made. However, denser disks develop 
a cool equatorial region and the disk densities and the resultant temperatures cannot be 
computed assuming a fixed, isothermal density structure. This result is confirmed by similar 
calculations for the disks surrounding the two later spectral types of Table [TJ However, rather 
than present similar plots and histograms for these models, we will proceed to the next section 
where a systematic comparison of the predicted density-averaged disk temperatures is made 
between fixed and hydrostatic models. 



To further investigate the reliability of the fixed density models, we examine how well 
they can predict the global, density-weighted, average disk temperature, defined as 



Here Moisk is the total mass of the disk, and the integral is over the total volume of the disk. 

Equation [12] was computed for 9 central stars, ranging in T e g from 13, 800 K to 30, 000 K 
in order to cover the full range of Be star spectral types. For each central star, 14 disk 



4.2. Density- Averaged Disk Temperatures 




(12) 
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Fig. 5. — The density-averaged temperature as a function of radial distance for six hydro- 
static models corresponding to values of p Q ranging from 5.0- 1(T 13 gem" 3 to 1.0 -lO -10 gem -3 
(see legend) for the BO model. 
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Fig. 6. — Comparison of the disk temperatures between a fixed (upper left) and hydrostatic 
(upper right) model with for p a = 1 • lCT 12 gcm _3 and n = 3.5 for the BO model. The fixed 
model assumes T Q = 13, 500 K. 
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Fig. 7. — Histograms of the density (left panels) and temperature (right panels) differences 
between fixed and hydrostatic models for p — 1 • 1CT 12 g cm" 3 and n = 3.5 for the BO model. 
The fixed models assume T Q = 9,500 K (top row), 13,500 K (middle row), and 17,500 K 
(bottom row). 
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models were computed for different choices of the base disk densit}|^|, p D , covering the range 
p = 2.5 ■ 1CT 13 g cm -3 to 5.0 • 10~ 10 g cm" 3 . All models assumed an n = 3.5 radial power-law 
index. In total, 126 hydrostatic disk models were computed. Figure [8] shows the results 
with the density- weighted, disk temperature expressed as a fraction of the stellar T eff . The 
(unweighted) average ratio at each T e s (shown in the figure by the filled squares) is well-fit 
by the quadratic relation 

/ j), " k * = 0.096 T 4 2 - 0.448 T 4 + 1.098 , (13) 



T. 



cff 



where T 4 is the stellar effective temperature in units of 10 4 K. Over the T e g range of the Be 
stars, the average ratio ranges between 0.55 and 0.65, suggesting some validity to the often- 
used approximation that a good estimat e for an (isoth e rmal) disk temperature is a constant 



fraction of the stellar T eS (for example, I Waters et al.l (119871 ) adopted < T Disk >= 0.8T cfT ). 



Similar ratios are found throughout the litera t ure, i ncluding some based on sophisticated 



modeling, such as that of ICarciofi fc Bjorkmanl (120061 ) who found a ratio of 0.6 for moderate 



density disks surrounding a B3 IV star. The current work finds that < Tbi S k > /T c s ~ 0.6 
is a reasonable fit to the average trend over the entire range of Be stars, with the quadratic 
fit of Eqn (TIB"]) representing a marginal improvement. However, the current work also makes 
clear that at each individual T e ff, the scatter about this average is large and depends in a 
systematic way on the disk density p a . Increasing the overall density drives the density- 
weighted temperatures to lower values due to the development of the cool equatorial zone 
noted in Figure [51 As shown in Figure ratios as low as 0.45 can be predicted for dense 
disks, and ratios as high as 0.7 for rarefied disks. 

To examine the accuracy of the fixed models in predicting the density-weighted disk 
temperature over a wide range of central effective temperatures and disk densities, the 126 
models were rerun but this time as fixed density models with T Q (see Equation [7] and discus- 
sion) chosen to be equal to the < T Disk > predicted by the corresponding hydrostatic model. 
Figure M compares the < T Disk > predicted by the fixed models to the hydrostatic results. As 
can be seen from the figure, the fixed models with an appropriately chosen T Q can reproduce 
the < TDisk > of the hydrostatic models to within ±500 K over a wide range of T e g- and 
p . This is particularly true for the models with low values for p Q , which scatter well within 
±250 K. However, larger deviations occur for the denser disks and, for T cff > 18, 000 K, the 
fixed models for the largest densities all predict hotter < T Disk > values by up to 1500 K. 

Despite this reasonable agreement, it is very important to next see if fixed models can 



2 The densities used were p a = 2.5, 5.0, 7.5 • 1CT 13 gem" 3 ; 1.0, 2.5, 5.0, 7.5 • 10~ 12 gem- 3 ; 1.0, 2.5, 5.0, 
7.5 • lO" 11 gem" 3 ; and 1.0, 2.5, 5.0 • lO" 10 gem- 3 . 
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Fig. 8. — Density- weighted disk temperature (Equation [T2l as a function of the central 
star's T eff . At each T eff , the open circles represent disks with different base densities, p a , 
with the size of the symbol signifying the size of p . The smallest circles correspond to 
p = 2.5 ■ 1CT 13 gcm~ 3 and the largest circles, to p Q = 5.0 ■ 1CT 10 gem" 3 . The filled squares 
are the (unweighted) average temperatures, and the solid line is the quadratic fit given by 
Equation $T3]). 
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Fig. 9. — Difference between the density- weighted disk temperature predicted by the fixed 
and hydrostatic models. For the fixed models, the T Q parameter was chosen to be the 
hydrostatic prediction for < T Disk >. As in Figure jS]), the open symbols represent the 
different base densities chosen for each disk. 
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be used to successfully predict more readily observed diagnostics such as Ha line profiles 
and infrared excesses. 



4.3. Observational Diagnostics 

To demonstrate more direct observational consequences of consistent vertical hydrostatic 
equilibrium, we considered a disk with a density of p Q = 5.0 ■ 10~ n gem -3 and n = 3.5. As 
noted in Section 14.11 this density is large enough that a cool equatorial region forms close 
to the star where the vertical temperature distribution is far from isothermal. Radiative 
equilibrium solutions were found for the fixed density structure of Eq. [TJ for several values 
of the T Q parameter. Observable quantities were computed for each of these fixed models 
to see if any could reproduce the observational predictions of the hydrostatic disk with the 
same p Q and n. This process was repeated for each of the three stars given in Table [1] which 
correspond spectral types BO, B2, and B5. 

The first observational diagnostic considered is the predicted Ha line profile correspond- 
ing to each model. The Ha line profile was obtained by solving the transfer equation using the 
hydrogen level populations predicted by bedisk. In all cases, a viewing angle of i = 35 de- 
grees was chosen (where i = corresponds to a pole-on star) and the equatorial velocity was 
set to 325 kms" 1 . This results in a projected stellar rotation velocity of v sin i = 186 kms^ 1 . 
An equatorial velocity of 325 kms -1 is approximately 70% of the critical rotation velocities 
given the parameters of Table [H Pure Keplerian rotation, with zero outflow velocity, was 
assumed for the disk. 

To represent the stellar disk, an LTE, photospheric, Ha line profile was used correspond- 
ing to the T e fj and log(g) of Table [H The profile for each element of the stellar surface was 
shifted by its projected radial velocity, resulting in a rotationally broadened photospheric 
pr ofile. For the calculat i on of the Ha emissivity and opacity in the disk, the Stark profiles 



of iBarklem fc Piskunovl (120031 ) were used 



Figure [10] shows the resultant Ha line profiles predicted by the hydrostatic and fixed 
models for each of the three spectral types considered. The fourth panel shows the variation 
with T a of the total Ha equivalent width of the fixed models. The equivalent widths of the 
consistent hydrostatic models are also shown. The Ha equivalent width increases with T Q for 
the fixed models and the increase is ~ 30% for a factor of two increase in T . Typically, the 
fixed model with (nearly) the lowest value of T Q is most successful in matching the hydrostatic 
prediction. This again reflects the influence of the cool equatorial region that develops in 
the higher- density disks. The Ha emissivity is controlled mainly by the high temperature 



-21 - 



"sheaths" above and below the equatorial plane in the inner disk. The lower T a values better 
represent the actual inner disk scale heights thus placing the hot sheaths closer the location 
predicted by the hydrostatic models (see Figure [3]). 

Another observational diagnostic is the infrared excess predicted by the models. Infrared 
excesses (expressed in magnitudes), relative to the underlying photospheric contribution, 
are shown in Figure [111 Again, the excess is found by solving the transfer equation along 
a series of rays threading the disk system. To represent the star, an LTE, line-blanketed 
stellar atmosphere was used corresponding to the parameters listed in Table [TJ Solutions 
were obtained at three wavelengths in the near- infrared, 2, 5, and 15 /im. The predicted IR 
excess for a series of fixed disk models with varying T a are shown in Figure [TlJ Also shown 
in the figure are the predictions of the hydrostatic models with the same p Q and n for each 
spectral type. For the earlier spectral types, the differences are generally small, within a few 
tenths of a magnitude, for plausible choices of T Q . However, larger differences are predicted 
for the latest spectral type considered, B5. Close inspection of this figure also shows that 
there is, in general, no single choice for T a that will reproduce the IR excess at the three 
wavelengths considered. Typically a higher T Q is required to match the excess at a shorter 
wavelength; this effect is particularly clear in the B5 model where the magnitude differences 
are largest. 



5. Conclusions 

In this work, we have compared predicted disk temperatures, densities, Ha line profiles 
and equivalent widths, and near-IR excesses between a set of Be disk models computed in 
consistent radiative and vertical hydrostatic equilibrium and a set of corresponding radiative 
equilibrium models that assumed a fixed density structure. Large differences between the 
predicted temperatures and densities can occur between the hydrostatic and fixed models 
when the density is large enough that the disk develops a cool equatorial region close to the 
star. In this case, there seems to be no choice for the single, isothermal temperature T that 
characterizes the fixed density structure that will yield a model matching the observational 
signatures of a consistent model. 

These conclusions are particularly important given the way such sets of Be disk models 
are used to compare to observations and extract disk parameters. Typically, a grid of models 
is computed for a wide range of p n with a fixed dens ity structure parametrized by a single 



temperature T Q (see for example ISigut &: Jonesl 120071 ). This grid is then used to compared 



to observations of, say, Ha line profiles or IR excesses to select the most appropriate disk 
parameters. As shown in this work, such a grid will contain a systematic error in the predic- 
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Fig. 10. — Predicted Ha line profiles for hydrostatic models and fixed models (assuming 
various T Q ) for three stellar spectral types (upper panels and lower left panel). The variation 
of the Ha equivalent width (in A) for these stars as a function of T Q for the fixed models 
is also shown in the lower right panel. The equivalent width predicted by the hydrostatic 
model associated with these stars is shown as an open circle. 
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Fig. 11. — Predicted IR excess for a series of fixed disk models computed by varying T Q for 
three stellar spectral types. Absolute IR excesses are shown on the panels on the left. The 
line styles and symbols correspond to wavelengths of 2 /im (solid line, circles), 5 /im (dashed 
line, squares), and 15 /im (dash-dot line, diamonds). The predictions of the hydrostatic 
models are shown as a dotted line with the wavelength in /im as indicated. The differences 
in the IR excess, fixed minus hydrostatic, are shown in the right-hand panel. The line styles 
have the same interpretation except that the dotted line in each shows the zero magnitude 
difference for reference. 
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tions for high p where a cool equatorial region in the disk develops; such models are poorly 
approximated by disks with a density structure fixed a priori. To extract the "true" distribu- 
tion of disk parameters, within the context of a thin disk in vertical hydrostatic equilibrium, 
a grid which consistently enforces both radiative and vertical hydrostatic equilibrium should 
be employed. 

We would like to thank Chris Tycner and Jon Bjorkman for many helpful discussions. 
This work is supported by the Canadian Natural Sciences and Engineering Research Council 
through Discovery Grants to TAAS and CEJ. 
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